<p> Something else </p>

d6cdc250-7249-415f-9c26-ee07fc1eb729 88 4 5 decimal

&[DATE] &[TIME] - &[FILENAME]

k gencm1 N 2 k matrix : i 1 k range N 1 i el i 1 - : N 2 i el k i - : 21 line for N 31 line :

Sanity checked 5 gencm1 012344321025 mat

M last M transpose 1 col rows :

nr nc v genm M nr nc matrix : i 1 nr range j 1 nc range M i j el v : for for M 31 line :

330 genm 00000000033 mat

M idx M transpose 1 col rows :

n k gencm n 1 ≡ k gencm1 M n 1 + 10 genm : i 0 k range N n 1 - k 1 + (i - gencm : w N idx 1 - : O 1 w 1 + i genm : N O N stack : M M N augment : 51 line for w M idx : M 1 n 1 + (2 w submatrix 41 line if 11 line :

without submatrix(...) n k gencm n 1 ≡ k gencm1 M n 1 + 10 genm : i 0 k range N n 1 - k 1 + (i - gencm : w N idx 1 - : O 1 w 1 + i genm : N O N stack : M M N augment : 51 line for w M idx 0 - : M 1 n 1 + (range 2 w range el 41 line if 11 line :

14 gencm 0123321024 mat

Only true for n=2 and any 'k'

24 gencm 11 mat 000001111222334012340123012010432103210210100315 mat 11 mat

n s weights d 2 s * 1 + : i 0 s range w i 1 + el 1 - (i^22 s * -^* d n + 2 i * - (d^i! (d n + i - (! (* / * : for w 31 line :

sanity 15 D 33 weights 0.043058311287478 0.090299823633157 - 0.054241071428571 0.008465608465608 - 41 mat

n s pts d 2 s * 1 + : i 0 s range M n s 1 + (i - gencm : M M 2 n 1 + 1 M idx submatrix : Pt i el 2 M * 1 + d n + 2 i * - / : 31 line for Pt 31 line :

Specific for n=1

n s 1 + (i - gencm

13 pts 0.8333 0.5 0.1667 13 mat 0.75 0.25 12 mat 0.5 11 mat 31 mat

n s pts d 2 s * 1 + : i 0 s range M n s i - gencm : M M 2 n 1 + 1 M idx submatrix : Pt i el 2 M * 1 + d n + 2 i * - / : 31 line for Pt 31 line :

Specific for n=2

n s i - gencm

24 pts 0.111 0.333 0.556 0.778 0.111 0.333 0.556 0.111 0.333 0.111 0.778 0.556 0.333 0.111 0.556 0.333 0.111 0.333 0.111 0.111 210 mat 0.143 0.429 0.714 0.143 0.429 0.143 0.714 0.429 0.143 0.429 0.143 0.143 26 mat 0.2 0.6 0.2 0.6 0.2 0.2 23 mat 0.333 0.333 21 mat 41 mat

n s pts d 2 s * 1 + : i 1 s range M n s 1 - (i - gencm : M M 2 n 1 + 1 M idx submatrix : Pt i el 2 M * 1 + d n + 2 i * - / : 31 line for Pt 31 line :

OBSERVE: 1. Nest 3 => kaput 2. Nests [1, 2] missing 1 col RHS

33 pts 0.125 0.125 0.125 0.375 0.375 0.625 0.125 0.125 0.375 0.125 0.375 0.625 0.125 0.375 0.125 0.125 0.375 0.125 0.625 0.375 0.125 0.375 0.125 0.125 0.375 0.125 0.125 39 mat 0.1667 0.1667 0.5 0.1667 0.5 0.1667 0.5 0.1667 0.1667 33 mat 0.25 0.25 0.25 0.75 0.25 0.25 0.75 0.25 0.25 0.75 0.25 0.25 34 mat 31 mat

<p> Something else </p>

split 0.955 0.864 0.773 0.682 0.591 0.5 0.409 0.318 0.227 0.136 0.045 111 mat 0.95 0.85 0.75 0.65 0.55 0.45 0.35 0.25 0.15 0.05 110 mat 0.944 0.833 0.722 0.611 0.5 0.389 0.278 0.167 0.056 19 mat 0.938 0.812 0.688 0.562 0.438 0.312 0.188 0.062 18 mat 0.929 0.786 0.643 0.5 0.357 0.214 0.071 17 mat 0.917 0.75 0.583 0.417 0.25 0.083 16 mat 0.9 0.7 0.5 0.3 0.1 15 mat 0.875 0.625 0.375 0.125 14 mat 0.833 0.5 0.167 13 mat 0.75 0.25 12 mat 0.5 11 mat 111 mat :

M Low High Scale High Low - (M M Min - M Max M Min - 106 -^+ / (* Low + :

j plot vx 1 split j row 1 el cols 2 range : vx vx 01 Scale : vy split j row 1 el transpose : XY vx vy augment : 41 sys :

10 plot 4 el 8 plot 4 el 1 plot 4 el 3 1 sys

n s pts d 2 s * 1 + : i 0 s range M n s i - gencm : M M 2 n 1 + 1 M idx submatrix : Pt i el 2 M * 1 + d n + 2 i * - / : 31 line for Pt 31 line :

Specific for n=2

n s i - gencm

Double integration algorithm f n s ∫ v 0 : Pt n s pts : w n s 1 - weights : i 1 w rows range P Pt i el : v v w i el P j col f j 1 P idx sum * + : 21 line for v 51 line :

t0 1 time :

x f e x 1 el^0.5 x 2 el (3^- (* :

accuracy doesn't increase wrt 's' ??? n 2 : s 7 : 12 mat

sanity ... 0.299449943475250 f n s ∫ 0.29944994347527

1 time t0 - 0.5 s *

=============================================

x y f e x^0.5 y 3^- (* :

t0 1 time :

acc 100 x y f y 01 x - int x 01 int 0.299449943465427

1 time t0 - 41 s *

sanity MCD wrt "kernel fnct" f(x) 0021 mat f 0121 mat f 1021 mat f 1121 mat f 14 mat 0.5 0.5 - 1.3591 1.3591 - 14 mat

n 2 : s 7 : 12 mat

part v 0 : Pt n s pts : w n s 1 - weights : i 1 w rows range P Pt i el : w P idx 31 line for P w P idx 13 mat 51 line :

part 0.3333 0.3333 21 mat 0.3634 0.8482 - 0.6768 0.2159 - 0.0247 0.0007 - 1.4898 106 -^* 71 mat 113 mat